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ABSTRACT 

We present the discovery of non-thermal radio and X-ray emission positionally coincident with the TeV 7- 
ray source HESS J 1 8 1 3—178. We demonstrate that the non-thermal radio emission is due to a young shell-type 
supernova remnant (SNR) G12. 8-0.0, and constrain its distance to be greater than 4 kpc. The X-ray emission 
is primarily non-thermal and is consistent with either an SNR shell or unidentified pulsar/pulsar wind nebula 
origin; pulsed emission is not detected in archival ASCA data. A simple synchrotron+inverse Compton model 
for the broadband emission assuming that all of the emission arises from the SNR shell implies maximum en- 
ergies of (30— 450)(B/10jU G)~° 5 TeV. Further observations are needed to confirm that the broadband emission 
has a common origin and better constrain the X-ray spectrum. 

Subject headings: acceleration of particles - supernova remnants - radio:ISM - Xrays:ISM - ISM: individual 
(G12.8-0.0,HESS J1813-178) 



1. INTRODUCTION 

Based on theoretical arguments, it has become widely ac- 
cepted that a significant fraction of Galactic cosmic rays with 
energies up to the "knee" in the cosmic ray spectrum at 
~ 10 15 eV are generate d in shell type supernova remnants 
(SNRs) (e.g. iBlandford & Eichlerffl987l) . While the detec- 
tion of non-thermal X-rays from the shells of several SNRs 
provides evidence that SNR sho cks are efficient accelerators 
of < 10 13 eV electrons (see e.g. Reynolds! IT996). direct evi- 
dence is still lacking that (1) SNRs can accelerate particles all 
the way to 10 15 eV, and (2) SNRs efficiently ac celerate pro- 
tons and heavy nuclei as well as electrons (see Pohll l200lt 
for a review). Powerful constraints on shock acceleration in 
SNRs can come from the detection of TeV 7-rays. However, 
these high energy photons can be produced by a variety of 
mechanisms, including inverse Compton (IC) emission from 
energetic electrons, and pion decay produced by the collision 
of energetic protons with surrounding dense gas. Addition- 
ally, pulsars, pulsar wind nebulae (PWNe), microquasars, and 
massive stars, have also been suggested as TeV counterparts, 
so it is unclear whether SNRs are in fact the primary source 
of cosmic rays < 10 eV. 

Currently only two convincing shell SNR/TeV as- 
sociati ons are known (SNRs G347.3- 0.5 and G266.2- 
1.2; lAharonian et alJ 120041 l2005al) . and observa- 
tions do not yet unambiguously distinguis h between 
hadronic and leptonic acceleration tLazendic et alJ l2004t 
Malkov. Diamond. & Sag deevl 120051) . However, with the 
advent of new sensitive Cerenkov imaging telescopes, there 
are brig ht prospects for discover ing additional TeV sources. 
Indeed, lAharonian et alJ (2005b) recently reported the dis- 
covery of eight new TeV 7-ray sources in the inner Galactic 
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plane (£ = ±30°) using the High Energy Stereoscopic System 
(HES S). Of the eight n ew TeV sources, six are listed in 
Aha ronian etafl (|2005b ) as being in close proximity to either 
an SNR, pulsar, or EGRET source (in some cases all three). 
However, these authors were unable to identify any plausible 
counterparts to two of the new TeV sources and suggest that 
they may represent a new class of "dark" nucleonic particle 
accelerators. 

In this Letter we present evidence that one of the two "dark" 
TeV sources, HESS J 1 8 1 3—178, is positionally coincident 
with a previously unidentified young radio and X-ray SNR, 
G12. 8-0.0. A simple model for the broadband emission is 
also presented. 

2. DATA AND RESULTS 

2.1. The Radio Source G '12. 82-0. 02 

The radio source G12.82-0.02 was discovered in a new low 
frequency Very Large Array ( VLA) 90cm survey incorporat- 
ing B+C+D configuration data (Brogan et al. 2005). The pre- 
viously unidentified radio source has a shell morphology with 
a diameter of ~ 2.5'. Analysis of archival VLA C+D config- 
uration 20cm data has confirmed this discovery, and indicates 
that the radio emission is non-thermal. G12.82-0.02 has inte- 
grated radio flux densities at 20 and 90 cm of 0.65 ±0.1 and 
1.2 ±0.08 Jy, respectively. Figure la shows the 90cm image 
while Figure lb shows a 3-color image of the region using 
90c m, 20cm, and 8 urn Sp itzer Space Telescope GLIMPSE 
data ( iBeniamin et alJ2003l) . It is notable that no distinct (non- 
diffuse) dust emission is coincident with G12.82-0.02, con- 
firming its non-thermal nature. HESS J 1 8 1 3—178 (positional 
uncertainty l'-2') is coincident with the radio shell (see Fig. 
la) and its size of ~ 3' is in excellent agreement with that of 
G12.82-0.02. No str ucture is evident in the p ublished 3' res- 
olution HESS image (Aharonian et al. 2005b). The radio and 
TeV sources are located ~ 8' from the Hll region complex 
W33 (see Fig. la,b, §3.1). 

Though not previously recognized as a distinct source, 
G12.82-0.02 is evident as a faint unresolved extension to 
W33 in the Nobeyama 3cm, Parkes 6cm, and Bon n 1 1cm sur- 
veys, with resolutions ranging from 3' to 4 . 3 ' llHavnes et alJ 
[19781 lHanda etltHlT^lReich et al][T990V The radio spec- 
trum of G12.82-0.02 (including the 20 and 90cm VLA 
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data, 4 m VLA B+C+D configuration data from Brogan et al. 
(2004), and the single dish data) is shown in Figure lc. 
The rms noise in the 20 and 90cm images are 5 and 2.5 
mJy beam" 1 , respectively, and the flux density uncertainty is 
(# independent beams) 5 x 3a. G12.82-0.02 is not detected 
at 4m; the flux density shown in Figure lc is a 5a upper limit 
(also see §3.1). To account for the extra uncertainty due to 
confusion with W33, 6a was used for the single dish uncer- 
tainties. The spectral index using these data (excluding the 4m 
non-detection) and a weighted least squares fit is -0.48 ±0.03 
(where S u oc v a ), similar to that of other small shell type 
SNRs (i.e. G 349.7+0.2 and W49B) without central pulsars 
JCreen 120041). 

Assuming a = -0.48, the 1 GHz surface brightness of 
G12.82-0.02 is Sjghz ~ 3.3 x 10" 2() W m" 2 Hz sr" 1 , and 
the radio luminosity from 10 7 — 10 1 1 Hz is ~ 3.6(d/d4) 2 x 
10 32 ergs s" 1 (d is the distance and d\ = 4 kpc; see §3.1). 
The radio morphology, radio spectrum, and lack of coincident 
mid-infrared emission (Fig. 1) provide convincing evidence 
that G12.82-0.02 is a previously unidentified shell type SNR; 
henceforth designated SNR G12.8-0.0. 

2.2. The X-ray Source AX J 181 3-1 78 

Analysis of a pointed 100 ks archival ASCA observa- 
tion toward this region (taken 1997 September) reveals a 
source, hereafter AX J 1 8 1 3—178, spatially coincident with 
both SNR G12.8-0.0 and HESS J1813-178 (see Figure 2). 
AX J1813-178 is detected by both the SIS and GIS ASCA in- 
struments, although the emission is contaminated in the GIS 
detectors by the bright X-ray binary GX 13+1 ~ 40' to the 
NE. Thus, we have confined our spectral and imaging analy- 
sis to the SIS data with a combined (screened) exposure time 
of ~ 64 ks. The X-ray source appears to be either unresolved 
or slightly extended (Fig 2). The X-ray peak nominally lies 
interior to the radio shell, but since the ASCA poin ting uncer- 
tainty can be as large as 1' (Gott helf etail 120001) . the ASCA 
image does not distinguish between X-rays originating from 
the center (i.e. a compact object) or the shell of the SNR. 

The X-ray spectrum is quite hard, with emission extend- 
ing to 10 keV and a strong cut-off below 2 keV; no lines are 
evident (inset Fig. 2). Thus it is clear that the emission is 
predominately non-thermal. The background-corrected SIS 
data (extracted from a 3' source region) consisting of « 2000 
counts in the energy range 0.5-10 keV were jointly fit to sev- 
eral absorbed spectral models including a power-law, and 2- 
component combinations of a power-law plus a (1) Raymond- 
Smith (RS) thermal plasma (with and without the abundances 
frozen at solar), (2) blackbody, and (3) Bremsstrahlung com- 
ponent. The data quality are such that all of these models give 
reasonable fits (reduced \ 2 ~ !)■ 

Statistically, a power-law plus RS model gives the best 
fit (RS kT ~ 0.22 ±0.1 keV) but the improvement is only 
4cr compared to a power-law alone; the emission measure 
of the RS component cannot be constrained. This result 
suggests that higher quality X-ray data are needed to deter- 
mine t he thermal contribution to the X-ray emission (see also 
Ubert ini et alJl2005T) . None of the two-component fits yield a 
significantly different absorbing column or photon index com- 
pared to the power law fit alone: Ng = (10. 8+23 x 10 22 cm" 2 
and T = 1.83~o42. The unabsorbed flux and luminosity of 
AX J1813-178 from 2-10 keV are - 7.0 x 10" 12 ergs cm" 2 s" 1 
and \.l{d/dn) 2 x 10 34 ergs s , respectively. These param- 
eters are consistent with AX J 1 8 1 3—178 being either a pul- 



sar/PWN or one of a small number of primarily non-thermal 
X-ray SNRs such as G266.2-1.2 and G347.3-0.5. 

3. DISCUSSION 

3.1. Distance Constraints 

As mentioned previously, SNR G12. 8-0.0 lies near the 
line of sight of the W33 Hll region/star formation com- 
plex. The radio recombination lines, molecular lines (in- 
cluding masers), and H I gas associated with W33 arise from 
LSR veloci t y co mponents at ~ +30 and ~ +50 km s" 1 (e.g. 
ISatolU977t Bieging etaD Il978t iGardnereTall Tl983). No 
Hi absorpti on features are evident tow ard W33 beyond ~ 
+55 km s" 1 ( Radhakrishnan et al. 1972). Since the tangent 
point velocity is ~ +170 km s" 1 in t his direction (as sumes 
Galactic center distance 8.5 kpc and iFich et alJfl989l rota- 
tion curve), W33 must lie at about the +55 km s" 1 near dis- 
tance of ~ 4.3 kpc. Similar to W33, H I absorption spec- 
tra from_lhe_Sojrfherr^ Survey (resolution 1', 
see lMcClure^Gri ffiths et al.l2005l) toward G12.8-0.0 show ab- 
sorption to, but not beyond ~ +55 km s" 1 , although the signal- 
to-noise is low. Thus, the distance to G12. 8-0.0 is likely to be 
> 4 kpc. 

Based on the derived radio spectral index (-0.48) the 
expected 4m flux density is ~ 2.7 Jy, suggesting that 
SNR G12. 8-0.0 should be readily detectable, yet it is not (Fig. 
lc, §2.1). This discrepancy could be due to the close proxim- 
ity of W33 if G12. 8-0.0 lies immersed in or behind the H II 
region due to free-free absorption (see e.g Kassim 1989). In 
interferometric low frequency images (v < 100 MHz) H II re- 
gions can appear in absorption against the resolved out diffuse 
Galactic plane synchrotron emission due to free-free absorp- 
tion. W33 is coincident with such a 4m absorption 'hole'; its 
morphology is similar to the diffuse 90cm and 8 fim emission 
(Fig. lb) and extends as far west as G12. 8-0.0. Thus, free- 
free absorption is the most natural explanation for the non- 
detection of G12. 8-0.0 at 4m, and confirms the H I result that 
the SNR must lie at the distance of or behind W33. 

The high column density derived from the ASCA data, Ng ~ 
10 23 cm -2 , also suggests that SNR G12.8-0.0 lies behind 
W33. Along nearby lines of sight, the integrated Gala ctic H I 
colum n density is only ~2 x 10 22 cm" 2 JPickev & Lockmanl 
1990), suggesting that a significant discrete source of ab- 
sorptio n must be in the foreground to the SNR. CO data 
( Dame et al] 120011) indicate that dense gas in the distance 
range 0-4 kpc (much of which is associated with W33) con- 
tributes a hydrogen column density ~ 8 x 10 22 cm" 2 , and can 
thus account for much of the X-ray absorption. This further 
argues for a distance of 4 kpc or greater for SNR G12. 8-0.0; 
we adopt a fiducial minimum distance of 4 kpc (dt) for the 
remaining discussion. 

3.2. The Probable Youth of SNR G12.8-0.0 

The small angular size of SNR G12.8-0.0 (2.5') suggests 
that it is quite young; certainly the high extinction toward this 
region would have prevented its detection optically. The ra- 
dius of SNR G12.8-0.0 is \.5{d/d A ) pc. If the SNR is still 
freely expanding at a speed v s its age is ~ 2%5{d/d^) (v s /vs)~ l 
years (v$ = 5,000 km s" 1 ). Shock speeds of v s = 5,000 
km s" 1 ar e expected from SNRs with X -ray synchrotron emis- 
sion (e.g. Aha ronian & Atovanlll999l) . If the distance is 4— 
10 kpc, use of the free-expansion age is justified since the 
swept up mass would still be a small fraction of the ejected 
mass. At greater distances the SNR is likely to have entered 
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the Sedov-Taylor phase. For example, at 20 kpc, the swept 
up mass is ~ 50 M© (assuming an ambient density n a = 1 
cm -3 and 10% helium abundance) and the Sedov-Taylor age 
is 2,520(w o )°' 5 (£'/£ , 5i)~ a5 years (where E is the kinetic en- 
ergy and £51 = 10 51 ergs). Thus, independent of distance, 
if G12. 8-0.0 is expanding into a typical density medium it 
is very young with representative age estimates in the range 
285-2,500 years. Excluding G12. 8-0.0, only six shell type 
Galactic SNRs with diameters < 3' are currently known, un- 
doubtedly due to observational selection effects since many 
more are expected (Greenl 120041) . Our discovery highlights 
the effectiveness of high resolution and sensitivity low radio 
frequency surveys in finding young Galactic SNRs. 

3.3. X-ray and TeV Emission from a Young Pulsar? 

Although there is excellent positional agreement between 
the radio, X-ray, and TeV emission (Fig. 2), it is currently 
unclear whether all these components originate from the SNR 
shell or if there is an as yet unidentified associated pulsar or 
PWN. No central radio nebula or point source is detected at 
our current spatial resolution and sensitivit y, and there are no 
know n radio pulsars within 16' of the SNR (Manchester et al. 
120051) . This field has been searched for pulsars as part of the 
Parkes multi-beam pulsa r survey, to a li miting 1 .4 GHz flux 
density of 0.2 mJy ( Manchester et alJ2001l) . corresponding to 
a pseudo-luminosity of ~ 3>{d/dn) 2 mJy kpc 2 . Recent deep 
pulsar searches toward SNRs have been able to either detect 
or place upper limits of ~ 1 mJy kp c 2 on the 1 .4 GHz p seudo- 
luminosity of central pulsars (e.g. Cam ilo et al.l l2002). Thus, 
the Parkes multi-beam survey is reasonably constraining, but 
a deeper search is needed. 

The X-ray luminosity and photon index inferred in §2.2 for 
AX J 1 8 1 3—178 are typical of young pulsars and their associ- 
ated PWN (e.g.[Possenti et al. 2002). In order to search forX- 
ray pulsations from AX J 1 8 1 3-1 78 we conducted aZ 2 search 
fsee iBuccheri et aLl fl983) on the barycenter-corrected ASCA 
high-bit mode GIS data (with 20 ks of data). No pulsed signal 
was detected using 585 counts between 4-8 keV. For peri- 
ods between 125 ms and 1000 s, we find an upper limit on 
the pulsed fraction of 38% (19%) for a sinusoidal (sharply 
peaked) pulse profile. Since most young pulsars rotate with 
periods below this range, and the pulsed signal may be a small 
fraction (< 10%) of the total pulsar plus PWN X-ray flux, the 
lack of pulsations is not very constraining. Deeper X-ray and 
radio observations of this source are needed to confirm or rule 
out a pulsar/PWN interpretation for AX Jl 8 13-178. 

3.4. X-ray and TeV Emission from the SNR Shell? 

The radio (SNR shell) and X-ray (origin uncertain) data fit 
smoothly onto a single power-law+rolloff spectrum, unlike 
the typical behavior of PWN, suggesting they have a com- 
mon origin (Figure 3). Assuming that the radio, X-ray, and 
TeV emission originate from the SNR shell, we have modeled 
the flux expected from synchrotron+IC mechanisms using a 



power law energy distribution modified by an exponential cut- 
off dN/dE w E'^expl-iE /E max )] (where a is the index of 
the electron distribution, a nd E mm is the maxim um energy 
of accelerated particles, see|Eazendic et al. 2004). Since the 
X-ray absorption column (and hence X-ray spectral index) is 
not well constrained by the ASCA data, two variations of this 
model are shown in Figure 3, one using the best fit X-ray Nu 
(red) and the other using the la lower limit to Nu (blue). To- 
gether these models encompass the likely range of parameter 
space. We have assumed a 15% filling factor for the magnetic 
field in the IC emitting region in order to obtain the best agree- 
ment with the HESS data for the 'blue' model; larger values 
yield lower IC components. The fitted broad-band electron 
spectral indices are a = 2.3 (red) and a = 2.0 (blue) - the for- 
mer is consistent with the best-fit photon index for the ASCA 
spectrum, while the latter is more consistent with the radio 
spectral index. 

The roll-off frequencies (v ro ) derived from the synchrotron 
spectra are 3.3 x 10 19 Hz (red) and 1.3 x 10 17 Hz (blue). 
The maximum energy is related to v TO and the magnetic field 
strength B by v ro - 1.6 x 10 16 (fi/10 /U G) (E max /\Q TeV) 2 
Hz, so that E max = 450 (B/10^ G)" 5 TeV (red) and E max = 
30 (S/lO/LtG)" - 5 TeV (blue). Thus, unless the magnetic 
field stren gth is much stronger than 10 fiG (see for ex- 
ample fVolk et al j |2005l). E mnr for this sourc e is potentially 
rather high (e.g. Reynolds & Keohane 1999). Additionally, 
for much of the parameter space range, the IC emission from 
CMB scattering alone is insufficient to account for the HESS 
data; the discrepancy is larger for larger values of B. It is 
possible a more sophisticated IC model incorporating scat- 
tering off starlight ( Aharonian et al J| 2005d) would bette r ac- 
commodate the HESS data. Indeed lHelfand et allJ2005l) find 
this to be the case assuming the SNR lies at the distance of 
W33, though currently it is only constrained to be at or greater 
than this distance. It is also possible that the TeV emission 
arises from the acceleration of nuclei and subsequent pion de- 
cay instead of IC if the SN R is embedded in a su fficiently 
dense ambient medium (e.g. lUchivama et al J 120051) . Detec- 
tion of resolved X-ray and TeV emission from the SNR shell 
would strengthen our suggestions that the radio, X-ray, and 
TeV emission have the same origin. Higher sensitivity X-ray 
data will also be crucial to better constrain the properties of 
the broadband emission. 
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FIG. 1.— (a) VL A 90cm image of G12.8 2-0.02 with 32" X 19" resolution. The cross indicates the position of HESS J1813-178 which has a positional 
uncertainty of 1' - 2' 1 Aharonian et al . 2005b). (b) Three color image of the same region as (a), with red=Spitzer 8 fim, green=VLA 20cm, and blue=VLA 90cm. 
Thermal emission is white to cyan or red in color, while non-thermal emission is darker blue, (c) Radio spectrum for SNR G12. 8-0.0. Solid hexagon symbols are 
from the current work at 20 and 90cm. The open hexagon symbol indicates the 5a 4m upper limit. The A symbols indicate flux densities from the 1 1cm Bonn, 
6cm Parkes, and 3cm Nobeyama single dish surveys i Reich et al. 1990; Haynes et al. 1978; Handa et al. 1987). The line (solid and dashed) shows a power-law 
weighted least squares fit to the radio data for v > 330 MHz (90cm). 
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FIG. 2.— ASCA SIS greyscale image of AX J1813-178 in the energy range 0.5-10 keV smoothed to 30" with 90cm contours from Fig. la at 28, 33, 38, 
43, and 48 mjy beam -1 overlaid (contour intervals are 2a). The black circle indicates the position and angular extent (3') of HESS Jl 8 13— 178. Inset ASCA 
SIS spectrum of AX J1813-178 extracted from 3' region encompassing the radio shell. The cross symbols show the data while the solid line shows the best fit 
power-law (dotted line) plus RS (dashed) model. The residuals in the form of the uncertainty a are shown in the bottom panel. 
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FIG. 3. — Fits to the broadband emission assuming that all the flux originates from the shell of SNR G12. 8-0.0. BPN stands for Bonn, Parkes, and Nobeyama. 
The diagonal black lines indicate both the uncertainty in the HESS flux measurements, and the fact that no spectral information has yet been published for the 
TeV emission. The two models indicated by the red and blue lines show the range of parameter space that best fit the data: the red model uses the spectral index 
from the best fit to the ASCA data Nn of 10.8 X 10 22 cm -2 (black X-ray spectrum), while the blue model uses the spectral index implied by the la lower limit 
toAf H of8.9 x 10 22 cm- 2 ( green X-ray spectrum). Both models include contributions from synchrotron (solid lines) and IC (dashed lines) mechanisms; we have 
assumed that the filling factor of the magnetic field in the IC emitting region is 15%. 



